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ABSTRACT 

K ■ We present broadband (radio, optical, X-ray and GeV) fits to the after- 

i-£h ■ glow light curves and spectra of three long-duration gamma-ray bursts (GRBs 

080916C, 090902B, and 090926A) detected by the Gamma-Ray Burst Monitor 
(GBM) and Large Area Telescope (LAT) instruments on the Fermi satellite. Us- 
ing the observed broadband data, we study the origin of the high energy emission, 
and suggest that the early-time GeV emission and the late-time radio, optical, 
C<") ■ and X-ray afterglows can be understood as being due to synchrotron emission 

^ ■ from an external forward shock caused by structured ejecta propagating in a 

wind bubble jumping to a homogeneous density medium. If the ceasing time 
^ ■ for majority of the energy injection is assumed to be close to the deceleration 

i ■ time of the forward shock, the structured ejecta with continuous energy injec- 

tion to the forward shock can well explain the early rising feature of the GeV 
7— i . mission from these burst, and the density-jump medium can account for some 

>• . certain plateaus or flares in the late afterglows. From our fits, we find that, on 

^ . one hand, the external shock origin of the GeV photons will make the optical 

depth have not significant contribution to the early LAT rising part, which will 
loosen strong constraint of lower limits of Lorentz factor. On the other hand, 
these Fermi-LAT events preferentially occur in a low-density circumburst envi- 
ronment, in which case the Klein-Nishina cutoff will significantly suppress the 
Self-Synchrotron Compton (SSC) radiation. Such an environment might result 
from superbubbles or low-metallicity progenitor stars (which have a low mass-loss 
rate at late times of stellar evolution) of type Ib/c supernovae. 



Subject headings: gamma-rays: bursts — gamma-rays: theory 
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1. Introduction 

Gamma-ray bursts (GRBs) are a kind of astrophysical phenomenon so luminous in the 
universe that their isotropic energies of 10 48 — 10 55 ergs are usually considered coming from 
extremely relativistic outflows with bulk Lorentz factors as high as 10 2 — 10 3 . The connection 
between long-duration GRBs and broad-li ne SNe Ib/c (low-metalli city) has been supported 
by some pieces of observational evidence (jWoosley fc Bloom! 120061 ) . The recently launched 
Fermi Gamma- Ray Space Telescope with the o n-board Gamma-ray Burst Monitor (GB M) 



and Large Area Telescope (LAT) instruments (lAtwood et al.l 120 09; Me egan et al.l 120091 ) in 
conjunction with the Swift narrow field instruments (IGehrels et al.ll2004j ) have opened a new 
era of understanding physical mechanisms of GeV photon emission in very energetic GRBs 
and their relation to lower-energy afterglow emission. Up to now, several novel yet tricky 
features have appeared during the whole period of observations and their complexities have 
challenged the anterior established theoretic models. 

The first feature is that the sub-MeV and GeV photons observed by GBM and LAT 
respectively behave distinctive spectral and temporal properties. The GBM light curves are 
nearly flat during the ma in episode of the pr ompt emission, then drop extremely rapidly [e.g. 
£-3.3 f or 080916C (lAbdo et al.ll2009ai )]. and eventually cease abru ptly. This can als o 

be validated from ~ 60% of all bursts detected by the Swift satellite ( lEvans et al.ll2009l ). 
For the LAT light-curves, however, a rise appearing in the early few seconds (< T 90 ) was 
followed by a single power-law decay without any cutoff up to hundreds of seconds after the 
GBM trigger, till below the monitor sensitivity, and decay of the late LAT emission wa s 
much shallower than the GBM-detected counterparts at the same times ( jZhang et al.ll2010l ). 
Furthermore, the LAT emission usually lags the GBM emission from a fraction of seconds 
to a few seconds. In addition, the spectral slopes of the GBM and LAT emissions are often 
different, e.g., the GBM data can be fitted with a Band spectrum that is composed of two 
smoothly-joining power laws, while the LAT data are often fitted by a power law with a slope 
intermediate between the two slopes of the GBM fit. These properties seem to indicate that 
the high-energy emission and low-energy emission detected by LAT and GBM have different 
origins. 

The temporal properties of the LAT emission have been studied, and some explanations 
differing in the source's dominant component have been proposed. A prevalent explanation 
is particle-dominated models. One inclination is the leptonic interpretation. Because of 
their distinctive light curve properties we mentioned above, it is highly probable that the 
MeV photons may be of an internal origin, while the GeV photons may be generated via 



(Kumar & Barniol Duran 


2009 


)• 


Wane et al. 


(2009, 


2010) 
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studied the Klein-Nishina effect on the high-energy afterglow emission and found that at early 
times such an effect strongly suppress the inverse Compton scattering of those electrons that 
produce the high-energy emission. Thus synchrotron radiation of t he electrons becomes a 
plausible mechanism. This conclusion was independently drawn by IZhang &: Pe'erl ( 120091 ). 
who suggested that the lack of a thermal component in the GBM spectrum of GRB 0809 16C 
is consistent with a relativistic Poynting-flux-dominated outflow. The explanation of syn- 
chrotron radiation of the electrons accelerated by a forward shock is fairly natural to account 
for both the observed delay of the > 100 MeV photons and long lasting of their emission. 
However, a thermal component deviating from the smooth Band spectrum function in some 
bursts (i.e. GRB090510, GRB090902B, GRB090926 A) seems to be bey ond a prediction of 
what is often invoked from the external shock model (IZhang et al.ll2010l ). Moreover, a rapid 
rise during the first few s econds (e.g. ^ t 6 of GRB080916C) is hard to be well explained 
within such a f ramework ( 



separate origin, 



Toma et al. 



bma et al.l 120091 ). Therefore, in addition to the hypothesis of a 



( 20091 ) assumed that GeV and MeV photons may come from the 
same region, but the onset of the high-energy may result from anisotropic inverse Compton 
scattering of an optical-thin expanding cocoon, delayed compared with the MeV emission. 
Nevertheless, in their calculation, this up-scattering cocoon is so short-lived that it could 
not account for the whole high energy emission. 



Another approach is the hadronic scenario. iRazzaque et al.l ( 120091 ) suggested that the 
MeV and GeV photons could be interpreted as the radiation of accelerated electrons and 
cosmic ray protons respectively, as well as the delay between them could come from different 
cooling time scales in a highly magnetized shock. In their framework, the model of GRB 
080916C is plausib le only when T < 500 and the jet opening angle ~ 1°. Alternatively, 
Asano et al.l ( 120091 ) considered another possibility of the photomeson cascade and proton 
synchrotron models, and provided their constraint in GRB090510. Notwithstanding, to 
reproduce the extra component around GeV with these models, the isotropic-equivalent 
proton injection luminosity is required to be larger than 10 55 ergs _1 . Such a large proton 
luminosity is a challenge for the hadronic models. 

All the works mentioned above inspire us to consider a plausible structured outflow, in 
which the bulk Lorentz factor of the initial shells tends to be lower than that of the late 
shells. This energy accumulation, therefore, would lead to an early rapidly rising light curve 
and transient soft to hard spectrum. 

The second feature of the three long GRBs with LAT data is th at some humps or 
flares appear in the light curves of the > 10 5 s low-energy afterglows dCenko et al.l 12010 



Swenson et al.ll2010l ). Neglecting this feature, iKumar fc Barniol Duran I (boioh fitted the 
multi-band light curves of these GRBs by assuming that a relativistic external shock sweeps 



-4- 



up an interstellar medium, and obtained reasonable physical parameters or parameter spaces 
accordingly. Although their fittings somehow favor the external shock model, humps (or 
sometimes flares) have indeed been observed on the late-time optical and X-ray afterglow 
light curves, which do not completely accord with a simple power-law but instead the decay 
after the humps is shallower than that before the humps. These observations call for a more 
meticulous consideration of the external shock model. 



It is noted that a density-jump medium proposed by lDai fc Lul (120021 ) has provided us 



2003; 


Tarn et al. 


2005 


; |jin et al. 


2009) 



stellar wind, subsequently encounters a density jump, and finally expands in a homogeneous 
medium. This interaction can produce an observed light-curve bump. 

In this paper, we show that the early-time GeV emissions together with the late-time 
radio, optical, and X-ray afterglows of GRBs 080916C, 090902B, and 090926A can be un- 
derstood as being due to synchrotron emission from an external forward shock caused by 
structured ejecta propagating in a wind bubble jumping to a homogeneous density medium. 
In Secf2]we include a set of observed broadband (LAT, XRT, UVOT) data on these three lu- 
minous long bursts. Then, our model is set up in Secj3j In particular, The structured ejecta 
can well explain the universal early rising feature of the GeV emission from these bursts, and 
the density-jump medium can account for some certain plateaus and flares in the late after- 
glows. Additionally, for the sake of verifying whether or not the lower energy (X-ray, optical, 
radio) emission originates from the same source as the higher energy (>100 MeV) emission, 
we discuss the effect of synchrotron self-Compton (synchrotron self-absorption) on the high 
energy (radio) emission, which, in the constrained parameter space estimated analytically 
from XRT and UVOT light curves, is proved in SecJUto have a small contribution to the flux 
density during the observed period. In this section, we find a reasonable set of parameters 
valid for most of the late afterglows. Our conclusions concerning a plausible central engine 
and ambient environment of bright, long GRBs are discussed in Secj5l 

Throughout this work, we adopt the convenience Q x = Q/10 x in units of cgs. 



2. Observations 



Among the 19 observed Fermi-LAT GRBs during the first 2 years' operation, GRB080916C, 
GRB090902B and GRB090 926A are 3 typical brightest long GRBs with abundant of spectra 
information ( IGranot II2010I ). 
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2.1. GRB080916C 



This burst, located at redshift z = 4.35 ±0.15 (IGreiner et al.ll2009[ ). is the first GRB de- 
tected by Fermi with high significance of photons at energies > 0.1 GeV. The isotropic energ y 
emitted from prompt emission can be esti mated as E^^ n = 8. 8 x 10 54 erg ( lAbdo et al.ll2009al ). 
At 00:12:45.613542 on 16 September 2008 f kbdo et al.ll2009at ). this GRB was triggered on by 
GBM with the duration T90 = 66s. Before To + 6s, the LAT light-curve shows an extremely 
steep rise F v j,a t oc t 6±0,5 , followed by a simple power law decay F u ^lat oc i- 1 - 33 ^ - 08 ^- 1 - 1 ^ - 1 
until ~ 1400s ( jZhang et al.ll2010l ). X-ray and optical photons were detected by XRT and 
UVOT since 17:11:28 16 September 2008 (T + 61ks). A steep decay (to ~ T + lOlks) 
continued with a plateau (to ~ T + 204ks) goes to a slightly shallower decay without break 
until 1.3Ms from the trigger T . During the period of T + 61ks— T + 1306ks, the sim- 
ple power law light curve and spectrum evolution show the flux density at X-ray band, 



oc t~ 



1.29±0.09,.-0.50±0.16 



V 



as well as for the optical band, F l 



v,opt 



^ £-1.40=1:0.05^-0.38=1=0.20 



(IGreiner et all 120091 ) . Besides, AGILE, RHESSO, INTEGRAL, Konus-Wind^and MES- 
SENGER all provide a plentiful of data information for the late afterglow ( IPerri et al.ll2008l ). 



Many groups hav e studied this GR B with the external shock model (IKumar fc Barniol Duran 



20101 : iGao et al.ll2009l ; IZou et al.l 120091 ). and suggested that the late time afterglow data can 
be used to extrapolate the early LAT data as well as be predicted from it, under the prefer- 
ence of circumburst density stratification. 



2.2. GRB090902B 



This burst is an exceptional case with redshift z = 1.822, whose speciality lied in 
both phases of prompt emission and early afterglow (lAbdo et al.ll2009bl ). Unlike majority 
of the other GRB events, excess emission exhibited in both low (< lOOkeV) and high (> 
lOMeV) band during the prompt emission. Moreover, the soft-hard-soft spectral evolution 
indicates two componen ts: Band function peaking at ~ 700keV + simple power-law with 
photon index T ~ 1.85 (jZhang et al.ll2010[ ). The isotropic energy release E lyiso = (3.83 ± 
0.05) x 10 54 erg ( ICenko et al.l 120101 ) . There are some explanations for origin of the high- 
energy power-law emission component observed in the LAT energy range (which accounts 
for ~ 24% of the total l OkeV to lOGeV fluenc e), including the hadronic origin [ either proton 
synch rotron radiation ( iRazzaque et al.l 120091 ) or pho tohadronic intera ctions (lAsano et al. 
2009|)] , or thermal emission from the jet photosphere (IRyde et al.ll2010l ). 



T his burst was triggered and loca ted at 11:05:08.31 on September 2009, b y the Fermi- 
GBM feissaldi fc Connaughton] l2009h and Suzaku-WAM flTerada et all l2009h with multi- 
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peaked duration T90 ~ 21s. After a rapid rise until To + 7s, the LAT band light-curve 
decays up to T + Iks with F UjL at oc t~ L4±0 06 ; including an energetic photon detected 
as high a s 11.16ln\« GeV w ithin the prompt emission phase and another 33.4jj^ GeV at 



Tn + 82s flAbdo et al.ll2009bh . T he Swift XRT (IKennea fc Strattal 120091 ; Evans et al.l 12009b 



Stratta et al. 20091 ) and UVOT (jSwenson fc SiegeT 20091 ) began concurrently target of op- 
portunity of the field of GRB090902B's fading source from 23:36 on 2 September 2009 and 
copious of data have been in hand from T + 0.5ks. A steep decay (to ~ T + 116ks) simulta- 
neously ended wi th a slight rise in both bands. It is highly possible that jet break occurred 
at ~ Tq + 553ks (jCenko et al.ll2010l ). During this whole period, the X-ray spectrum evolved 



as F,, x oc t~ 



1.36^10.03^-0.^0.13 



, and for the optical band, F l 



u,opt 



^-o.sgio.os^-o.yeio.o? 



In ad- 



dition, VLA began to observe the afterglow since 3 September 2009 (jChandra fc Frailll2009l ) 
until 5 months later at 8.5GHz and 4.8GHz. 



Using the forward-reverse shock and constant density model, ICenko et al.l ( 120101 ) sug- 
gested the afterglow of this burst is better fitted in X-ray and radio bands than in the 
optical/UV band, except for the first point at ~ 10 4 s. Parameter constraints without taking 
LAT data into account suggest a low circumburst den sity and a large kine tic energy, which 
is possible in the narrow opening angle ~ 3.4°. Later, iLiu fc Wangl (120101 ) have provided a 
meticulous discussion about the two-component forward-reverse external shock in a mono- 
tonic circumburst environment and their calculation was well fit to the four observed bands, 
except for the very early rising part of the LAT light curve. 



2.3. GRB090926A 



This burst is the GRB detected at z = 2.1 062 by Fermi with photons as high as ~ 20 



GeV at 26s after the trigger (jUehara et al.ll2009l ). The prompt emiss ion time-integrated flux 
indicates its isotropic emission energy of E 1 ^ iso = 2.1^Q Qg x 10 54 erg (jZhang et al.ll2010l ). 



Since the trigger at 04:20:26:26.99 on 26 September 2009, GBM and LAT started oper- 
ation and found that multiple pulses with total duration T 90 ~ 20s can be well modeled by 



a Ban d function spe ctrum in the prompt em ission, along with the Suzaku-WAM dNoda et al. 



2009 ). Konus-Wind flGolenetskii et al.ll2009h . and RT-2 on CORONAS-PHOTON (IChakrabarti et al 
20091 ). The light curve rose slightly for quite a few seconds (to ~ T + 16s), and decayed 
F u ,lat £- 2 - 05±0 - 14 z/~ L26 -o.22. A f ading X-ray counte rpart was observed by XRT since 17:17 
on 26 September (T + 46.7ks) JVetere et al.l l2009h . then PROMPT, SMARTS detected 
photons in several optical filters but no radio (5.5GHz) source under the limit of 1.5mJy 
was detected up to 1 October 2009. As shown on the light curve, a plateau (T + 51.4ks 
~ T Q + 92ks) overlaps the simple power law decay that is slightly steeper than the fading 
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supposed to be extrapolated from the early trace. No break up to lOOks from the trigger 
T was observed. During the period from T + 46.7ks to T + 149ks, two variabilities, are 
suggested to be flares (the first, at ~ 70ks — 95ks with St/t ~ 0.35; the second is slight, 
at 195ks-260ks, with 5t/t ~ 0.28) overly on the simple power law ligh t curve in both X-ra y 



band (F v> x oc t 



-1.40±0.05 



v 



=+0.3 
->-(). 2 



)and optical band (F Vt0pt oc t 



-1.01] 



(ISwenson et al.ll2010l ). 



Cenko et al.l (120101 ) and iRau et al.l (120101 ) studied the later afterglow (X-ray, optical) 



of this burst, provided a certain parameter space, confirmed the second flare, and indicated 
that the jet break occurs around ~ 2 Id. However, due to lack of radio data, the constraints 
cannot be narrowed down. Meanwhile, as for a commo n rebrightening plateau in 5 optical 



bands, they sugg ested a density-jump circumstance (e.g.. lDai &: Lull2002l ; lLazzati et al 



2002 



Tarn et al.l 120051) or a smooth injection of energy in to the forward shock from the central 



engine (e.g. 
the early higher energy emission. 



Dai &: Lulll998al ; iRees fc Meszaroslll998l ). and called for a detailed analysis for 



Model 



The behavior of the GRB afterglows shed light on the external shock model ( IPanaitescu fc Kumar 



20021 ). A collimated (6 < 10°), ultrarelativistic outflow of matter and/or radiation is driven 
by a central engine. Some certain dissipative processes within the outflow give rise to the 
prompt gamma-ray emission, with a fraction <^ 7 = E ltiso / (E lyiso + EK,iso) of the total rela- 
tivistic energy converted to high-energy radiation. 

A self-consistent result of the multiband (X-ray, optical, radio) afterglow data coor- 
dinated with the prompt emission relics producing > lOOMeV photons, may support the 
external shock origin of LAT photons, which is the same as the later lower-band emission. 
Meanwhile, due to the difference between the light curve slope of LAT and GBM, we take 
no consideration of the GBM emission in this piece. Moreover, the particularities of GRB 
090902B we mentioned above (soft-hard-soft spectral evolution, flares) may evoke doubts 
about the conventional external shock + monotonic circumstance model (we will discuss 
it in Secj5]), but enlighten our deliberate consideration on amendment to structured ejecta 
sweeping up the density-jump medium. 
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3.1. Dynamics 

3.1.1. Before the Deceleration Time 

Energy can be injected by a central engine continuously during a period (majority of 
which is within the first main pulse of prompt emission, before t/J3). We here suggest, within 
the framework of the collapsar model, an increase of the Lorentz factor of the ejecta in active 
time of the central engine may due to two reasons. First, angular momentum of the accreted 
fall-back matter would spin up the central compact object and thus its rotational energy loss 
could give rise to an increase of the ejecting luminosity. Second, earlier-ejected shells, when 
breaking through the stellar envelope, may suffer from more massive baryon contamination, 
and thus later-ejected shells may propagate in a tunnel (the later the cleaner). Therefore, 
the Lorentz factor of the ejecta head may increase with time. For simplicity, we assume the 
bulk Lorentz factor of the front materials distribute as a power-law function of time: r oc t x 
for t < tf p . If x > 0, the Lorentz factor of the shocked matter increases, being due to energy 
injection. 

We denote M e j as the accumulated mass ejected by the central engine at some certain 
time, and rj as the bulk Lorentz factor of a blast wave at this certain time, which is defined 
when the mass of the surrounding matter swept up by the blast wave, M sw , is equal to M e j/rj. 
After this time, the bulk Lorentz factor presents an obvious deceleration, contributing to a 
peak of the LAT light curve. This time is called the deceleration time td ec , which is assumed 
in this paper to be around the stop time of effective energy injection. This assumption 
is reasonable because energy injection to the forward shock is effected by the circumburst 
medium. Thus, both the apparent energy- injection time tj p and the deceleration time td ec 
should be shorter than the central engine ceasing time Tgo. 

The shock propagates a distance SR ~ 2(1 + z)~ 1 Y 2 cSt during a small observed time 
St. Assuming the proton nu mber density of a medium n oc R~~ k , thus for a homogeneous 



medium (k = 0), n = const (ISari et al.l 119981 ) or for a wind medium (k = 2), n = AR 



-2 



1 As Fig. 2, 10, 11 in IZhang et al.1 (j2010h shown, majority of the > 100 MeV photons in the prompt 
emission are released during the first largest pulse (at - 10s - 10s - 10s for GRB 080916C, GRB 090902B, 
and GRB 090926A, separately). Comparing this time to the peaking time of the LAT lightcurve (~ 6s, ~ 
7s ~ 16s for GRB 080916C, GRB 090902B, and GRB 090926A, separately), and the T 90 (~ 66s ~ 21s ~ 20s 
for GRB 080916C, GRB 090902B, and GRB 0909 26A, separately ) , we d efined t fp , before when energy has 



been largely injected. Of course, as mentioned in iMaxham et al' ( 2010h . whole of the energy injection last 



longer, while it will dwindle after the first largest pulse. 
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10 35 - 5 A 35 . 5 R~ 2 f lDai k Lulll998bl : IChevalier fc Lill200Q[ ). In both cases, we obtain 



_ gx+l)(l + *) 

''dec „ o -'l-dec- 



2l] 2 C 



The deceleration radius can be written as 



^ f 3 — k -Ek.isc 
60 \4:7iAm p rfc 2 



l/(3-fc) 



(2) 



Therefore, < t^ec) = Rdec(t/tdec) + ■ According to ISari I (119971 ). in order to check 
whether the the unshocked ejecta is assumed to be equal to that of the shocked matter, let 
A = ctf p (l + z)~ l represent the thickness of the shell, and / be the density ratio between 
the preshock fluid in the shell and in the circumburst surrounding (k = 2), given by 



./ 



1 + x)(l + z)E K<iso tf p 



H+l 



4nAR 2 - k m p c 3 r] 2 t"+ 2 



t<t fp 



(3) 



This may bring about a correction for the GeV rising part under different conditions, i.e., in 
the thick shell case, /<r 2 , a reverse shock is relativistic, the Lorentz factor of the shocked 
fluid should be changed and the contribution of the reverse shock to the flux density should 
be considered, while in the thin shell case, / ^ T 2 , this change is negligible. However, no 
matter how the thickness of the shell is, tf p < t^ c . In the following analytical calculations, 
we simplify that tf p ~ t dec (which can be observed from the figures of temporal flux density) 
and assume that the self-similar condition is established after the deceleration time. 

Before the deceleration time, energy is continuously injecting into the former shells. 
The mass ejected later catches up with the earlier ones, and V can be simplified as the bulk 
Lorentz factor of these materials combination. In the very beginning, it is possible that 
the bulk Lorentz factor is too low for majority of high-energ y gamma rays to escape from 
the MeV background. According to predominate hypothesis (IGranot Il2010[ ). only when the 
optical depth of the 7-ray absorption r 77 < 1, the bulk Lorentz factor i s large enough to 
radiate almost completely. In this piece, we use sketchy in lZou et al.l ( 1201 ll ) to check whether 
the absorption play a dominate role in certain parameter spaces. 



3.1.2. In the Density- Jump Surrounding 

Because of the association of long GRBs with star forming regions, it is highly possible 
that massive stars are still embedded in a cloud when giving birth to GRBs. Low density 
bubbles are created by stellar winds from GRB progenitors, whose sizes and densities strongly 
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depend on the initial ambient density. Therefore, a density jump occurs at the boundary 
between the wind bubble and the outer cloud, w hich may result in a bump/flare and shallow 
decay of the light curve during the later period (jTam et al.ll2005l ). 



From the moment tbreak at which the blast wave reaches the boundary, to the termination 
when it comes out into the interstellar medium at t e nd, the number density of the medium 
changes as 



3 x 10 35 A 35 . 5 R~ 2 , 



n 



AR~ 



3 x 10 35 [1 



log_f?.-log_B b 

1' 



re ah 

n aftcr = 3 x 10 35 eA 35 .5iCak> 



(f - 1)}A 35 . 5 R 



2 

break ' 



R < -Rbreak, 
-Rbreak < R < -Rend, 
R > -Rend, 



(4) 



where £ is the ratio of the number density in front of and behind the jump. 



It has been usually assumed t hat the afterglo w emission from the blast wave is nearly 
adiabatic ( jMeszaros fc Reea 119971 : ISari et al.lll998l ). so the total kinetic energy of the rela- 
tivists shock is constant. However, if we take radiative energy loss into consider ation, the 
radiat ion efficiency of the blast wave can be given with the electronic fraction e e by lWu et al. 
f)2005h 



-et 

fVmN(p-2)/2 



Vrr, < V„- 



(5) 



In such a situation, energy loss is significant in the fast cooling and gradually dwindle to the 
quasi-adiabatic case (i.e., the slow cooling phase). 

Hereafter, a denotation m — (3 — k)/(l — e) can analytically result in the bulk Lorentz 
factor evolution, 



77 ( 



'dc 



-m/(2m+2) 



Consequently, the Doppler factor becomes T>- 
emitting latitude angle). 



1 



t ^ ^dec, 
t ^ ^dec- 

2r 



r(l-/3cos6») 



i+r 2 e 2 



(6) 

2r (where 9 is the 



When the blast wave reaches the density boundary, the observed light curve will not 
behave as an abrupt 



ump but a smooth rise or pla teau instead, due to the curvature effect 
( iFenimore et al.l 119961 : iKumar fc Panaitescu I l2000h . Although this is trifle to be accom- 
plished in analytic modeling, a simplified power law distribution of medium density in the 
period of tbreak ~ tend can account for the same light curve feature (See Fig. Hla). Succeed- 
ingly, the blast wave radius can be written as 

t \2^+l 



R 



-Rdec(rM 

''dec 



(l-e)/(4-fc-e) 



t — tdecy 
t — ^-dec- 



(7) 



- 11 - 



3.2. Radiation 



In the comoving frame, with an electronic fraction e e and magnetic fraction e# of the 
post-shock thermal energy density, which may either vary or not with different media. De- 
note f p = 6(p — 2)/(p — 1), the magnetic field intensity B = Y(32u:e B nm p c 2 ) 1 / 2 . The 
cooling, minimum, and maximum electron Lorentz factors are 7c = 7syn)C /[l + ^(7c)] = 
67rm e c(l + z^arB^t)- 1 /[l + F( 7c )], 7m = (l/6)/ p e e (m p /m e )r, and 7max ~ 10 8 /\/S, 



while the synchrotron characteristic frequencies v c 



can be calculated by z/j 



/ -ffVeB/[2iTm e c(l + z)\ (where i = c,m, max), where F( 7 *) is the Compton factor for elec- 
trons with the Lorent z factor 7*, defined as the ratio of the self-inverse Compton (SSC) to 
synchrotron emission ( jNakar et al.ll2009l ). 



SSC 



(7*) 



p. 



</H (7* 



When Klein-Nishina (KN) effects are unimportant, i.e., the IC scattering of 7 * electrons with 
synchrotron photo ns are in the Thomson scattering regime, F( 7 *) can be simply derived by 
( Sari fc Esinll200ll ) as (a/1 + 4ee e /eB — l)/2, which is independent with 7 *. However, for 
high-energy electrons whose KN effect becomes import ant, Y('y^) depends on 7 * and can be 
expressed as analytical solutions in different cases from I Wang et al.l (]2010l ). 



The total number of shocked accelerated electrons N e = 4irnR 3 /(3 — k). The peak 
spectral power of a single electron is P Ujm _& x = (TTm e c 2 VB / (3e) , and the peak flux density 
at certain time of the afterglow is F u<raax = (1 + z)N e P u<raax / \AtiDl 2 ). Here, the luminos- 

z 

ity distance D L = (1 + z)H ^ 1 c j[Q m (l + z) 3 + Q k (l + z) 2 + Vt^ l / 2 dz (here we adopt a 



Cfpr 


standard ACDM co smology with H = Tlkms^Mpc- 1 , Vt m = 0.27, tt k = 0.0, fi A = 0.73 
JSpergel et alil2007K 



For a lower energy band (e.g., radio), synchrotron self-absorption may play a crucial 
role. The self-absorption optical depth r v , m or (r^ c ), according to the definitio n r v = 



becom es r u , 



c o(p-l) enR 



and r v , 



CO 



f4, where c « 10.4(p + 2)/(p + 2/3) flWu et al. 



(3-fc) B 7m 5 (3-fc)B 7c 5' 

20051 ) . Succeedingly, the self-absorption frequency v a can be written as 



r. 



T, 



3/5 

.1 v u 

2/(p+4) 



;/./ 



T 2/(p + 5W )1/(p+5 )^ 



< 1, 

1 < T VJL < (^) (p+4)/2 , 
2fc)(p+4)/2 <Tvt< (l 



(9) 



)(p+5)/2f£LU/2 
' v v h > 
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Then, the SSA power-law light curve yields ( ISari &: Esinll200ll ) 



(ii)-(9-l)/2 ) 
f^.V(9-l)/2fii 



-p/2 



V < U a , 

V < v h 
Vi<V< v h , 

V > v h . 



(10) 



Here v\ = min(z/„ 



max^ m , is, 



q 



2 when v c < is m , while v c > is m , q 



P- 



It is noted that before the optical depth becomes thin, the annihilation effect of high- 
energy photons can affect their flux, so the observed flux can be estimated as F ^ tU = e~ T ~ n F u . 
After T 90 , F ohjU = F v . 

For a high energy band, synchrotron self-Compton (SSC) radiation might contribute 
to the observed flux, the characteristic frequencies f/c,c, ^ic,m, ^/c*,max, -F/o,max can be 
calculated by vic,i = 2Uijf (where i = c, m, max), and Frc,u,- m^ = (l^/^5)aTRnF Utm3iX . Thus 
the SSC power-law light curve yields (IGupta fc Zhangl 120071 ) 



uJC 



/C,v,max 



- via,l ' 
v ' 



,1/3 



-{q-l)/2 

(^I£±\-(q-l)/2 { 
^ VIC,1 ' 



-p/2 



V < VlC,l, 

vie, i <v < vic,h, 

VlC,h < v 



(11) 

Here v IC ,i = mia(u IC , m , vic,c), vic,h = max(z/ /c - im , z/ 7CjC ). q = 2 when u IC , c < vic,m, while 
vic,c > vic,m, q = P- 

However, due to the KN effects on optically thin spectrums , Equations ([TU]) an d (fTTj) 
should be amended as different cases mentioned in equations of iNakar et al.l (120091 ) . This 
calls for the definition of electron KN Lorentz factor: = m e c 2 T/[hu syn ('yi)] (i = c, m). We 
have, if the synchrotron photons emitted by electrons with Lorentz factor larger than 7; [i.e. 
with frequency > v S yn{li)]-, they cannot be upscattered efficiently by electrons with Lorentz 
factor larger than 7, (which is above the KN limit). 

In the case of a jet, a simple assumption is that the jet decelerates without sideways 
expansion. After the jet break time tj et at which the jet angle 8j 
is assumed to be along the line of sight), the observed flux F U j Ct = 



^ 1/r (where the jet axis 
F u (t/t jet 



,(3-fc)/(4-fc) 



Comparing the flux density estimated analytically above with the observational multi- 
band values, certain physical parameters can be constrained. Whether the parameter space 
is reasonable or not can on some extent verify the feasibility of our external shock model. 
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4. Fitting Results 

4.1. Constraints on Parameters 

Based on the analytical model we mentioned above, our procedure of further narrowing 
down the physical parameter space in the self-similar phase includes the following steps: 

In the wind environment: 

(1) Accordi ng to the expression of F u oc t~ a z/ -/3 in different spectral regimes given by 



Wu et al.l (120051 ) and the observed high-energy, X-ray and optical temporal and spectral 



indices, (p,e) can be roughly estimated (see Tab. [TJ). 

(2) For simplicity, e ~ e e in the fast cooling regime and then e decays slowly (where 
we denote t cm as the time when v c = u m , and obtain that e decays as e ~ e e (t/icm) 2 ~ p in 
the slow cooling regime (i.e., setting j = —m/(2m + 2), we have, after t cm , F oc P , then 
v m oc t~ (j+1 ) and v c oc t^ +l \ and equation (JSJ) provides e oc e e (u m /u c Y P ^ 2 ^ 2 oc t^ 2 ~ p ^ j+1 \ as 
long as e < 2/3, j + 1 ~ 2/3). 

(3) Because of the definition for the deceleration time [set k = 2 in Equations ([T]) and 
()2])], r] can be expressed with x, -EK,iso and n (or A 35 5 ). This can estimate whether the shell 
thickness is thin enough to neglect the reverse shock's contribution to the total flux density. 

(4) Hereafter the free parameters are x, E K ^ iso , ee, and n (or A 35 5 ). The observed 
flux density in different wavebands (radio, optical, X-ray, >100 MeV) can be combined and 
eventually simultaneously be solved. In this paper, the temporal indices of the X-ray and 
optical light curves concurrently evolve according to equation fTlO|) and thus e# can therefore 
be expressed with x, -EK,iso and n (or A 35 . 5 ). 

(5) To justify whether the early high energy LAT photons come from the same region as 
the late afterglow, what we need to do is to replace es in the form of x, -EK.iso and n (or A) 
from step (4). To confirm whether SSC make a significant contribution to the GeV afterglow 
or Klein-Nishina effect changes the synchrotron and SSC spectra, we need to discuss the 



value of Y{pf^) in the slow cooling regime o r Y(~f m ) in the fa st cooling regime ( jNakar et al. 



20091 ). therefore Yi^Gev) can be estimated (jWang et al.ll2010l ). If no solution is found in the 



range of (10 53 ergs ■< E KMo ■< 10 56 ergs and 10" 2 ■< A 35 . 5 ^ 10 2 ), the >100 MeV photons 
may come from the other component, vise vesa. 

(6) From discussions about the parameter space of (A 35 5 and x) according to observa- 
tional statistics, one can obtain several best fittings. In addition, The ratio £ of the wind and 
homogeneous medium density, e e and e# at early and late times may have slight changes, 
which can also be legally adjusted in practise. 
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4.2. Best Fittings 



We put the above procedure into operation, and best fitting parameters are shown in 
Tab. [21 The evolution of character istic frequenc i es for each burst is expected as follows. 
Here our LAT data are taken from IZhang et all (120101 ) and X-ray data were reported by 
Swift fevans et all I2007L l2009h . 



From the lightcurves, we found that from tde. r . on, ao e y,x,opt > 0, PceV,x,o P t > 0. Ac- 
cording to the analytical expression from IWu et all (120051 ). before the blast wave reaches 
the density jump, GeV photons at early decaying time (e.g., t dec < t = tj < 10 3 s) has two 
possibilities: 



F S yn,GeV ^ 



3p-2-(p-2)e -f 
f 2(2-e) u GeV 

3p-2-(p-2)e , p-2 
t 2(2- e ) T" . ■ 



'V. 



GeV 



VGeV > V m (t) > U c (t), 



while the lower energy bands (X-ray, optical) at t — tjj > 10 4 s, may lie in one of the following 
slow cooling cases: 



3p-l-(p-l)e -(P- 1 ) 
2(2- £ ) Vx 2 



3p-l— (p— l)e -(p- 1 ) 

F opt xt SP=S) v opt 2 ,F x ccf 

3p— l-(p-l)e -(p-1) 3p-2-(p-2)e , p-2 -p 



F opt oc f 



3p-2-(p-2)e 1 p-2 ~P 



2(2-e) 



3p-2-(p-2)e , p-2 ~P 
2(2- £ ) + 4-pjy^2 ^ 



VX > V c {t) > V opt > V m (t), 
V X > "opt > Vc{t) > 



for 090902B radio band, at t — tm when j5 ra d%o > (i.e.i/ r<M jj > i/ m ), it is similar to the 
optical spectrum and lightcurve. While (3 ra dio < 0, 



Fradio OC t 3 < 2 ^^ ^ , 



radio- 



4.2.1. GRB080916C 



Our optical photometry of this burst is taken from iGreiner et all ( 120091 ). First of all, if 
the X-ray and optical afterglow shares the same source as the high-energy photons detected 
by LAT, their spectral indices f3 b )X = 0.50±0.16 and /3 ^ op t = 0.38±0.02 should be consistent 
with the GeV photons, a b,LAT = 1-33 ± 0.08 and Pob,LAT = 11 ± 0.1. This is preliminarily 
valid only in the case of vcev > max[i/ c (£j), v m {ti)\ and u c (tn) > v x > v opt > v m (tn), when 
2.0 < p < 2.32 and < e e < 0.48. 

Consequently, the peak of the LAT light curve lying at td ec ~ 6s can be considered as 
the deceleration time of this burst. Therefore, equations ([1]) and d2J) lead to the relationship 
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among the isotropic kinetic energy £k,55> the bulk Lorentz factor 773, and the wind parameter 
A 35 , 5 at the very moment: r] 3 ~ OAOE^A^il + 2x) 1/4 . 

Equation © provides, at t < t dec , f /T 2 = 72 exp(5.37x)(l + x)(l + 2x)- 1 r 3x ~ 2 > 1. 
Thus the reverse shock is nearly Newtonian. 

Subsequently, when the blast wave is embedded in the wind bubble, we assume that the 
typical value e{t{) ~ 1/3, while e(tjj) ~ 0.0, and the characteristic frequencies evolve as 

u c (tj) = 1.21 x 10 9 Hz4^ 5 A 3 - 5 2 5 e B 3 /2 (l + 2x)- 5 /V/ 5 [l + Y( lc ,tj)]- 2 , (12) 

i/ c (* 7J ) = 1.03 x 10 9 Hz4^ 5 A 3 - 5 2 5 es 3 /2 (l + 2x)- 5 / 2 t 7 // 2 [l + F( 7c , t 7/ )]" 2 , (13) 

v m (ti) = 8-66 x 10 23 Hz^e 2 ^(l + 2x) 3 / 2 t 7 ~ 8 / 5 / P 2 , (14) 

v m {t u ) = 7.39 x 10 23 Hz4( 2 5 e 2 e^ 2 (l + 2x) 3 / 2 t 7 r 3 / 2 / p 2 . (15) 



Setting p = 2.3 and e 6j o = 1/3 as a trail, we only consider the synchrotron emissivity, 

iV,Gev(13s) = 18.5/xJyEgf e%£(l + 2x) 9 / 40 [l + Y( lc , 13s)]" 1 . (16) 

We found that equation (fTBj) is approximately the same under both fast cooling and slow 
cooling condition, so do with the following equation ( 124)) and (151)1 . 



*W(1.01 x 10 & s) = 1.76 x lOVy^f ^W^H 1 + 2 ^ 



33/40, 



^59/40 



7.73 x 10">Jy. (17) 



From equation ([171), we have e B ,o = 1-96 x 10~ 8 E^ 55 A 3 ^ 33 (1 + 2x)~ 9 / 5 , so t cm < 13s in the 
case that 10~ 2 ^ ^35.5- Therefore, u m (13s) < z/ c (13s) and F S2/nj G e v(13s) < 0.83/iJy, which 
yields 



5.83E Kj55 Ai /10 (l + 2X) 1 / 10 < 1 + F( 7c , 13s) < 54Ek,55^35 5 5 (1 + 2*0 



-1/5 



From lWang et al.l (120 lOl ). we have, in the wind, if tj > 13s, Y(^ c , i 7 ) < F(7 C , 13s). This have 
been shown in Fig. Hlc. Hence, as long as A 35 . 5 (l + 2x) _1 t 7 < 20.2 (e.g.l0~ 2 ^ ^35.5(1 + 
2x)~ x < 1.56 at 13s or 10~ 2 < A 35 5 (1 + 2x)~ x ~ 0.02 at 10 3 s), 



= 3.0x 10- 3 E K 1 55 A 2 / 5 5 t 7 3 / 5 (l + 2x)- 2 / 5 [l + F( 7c ,t 7 )] < 0.18^ 5 5 t 7 3 / 5 (l + 2x)- 3 / 5 < 1, 

le\tl) 

(19) 

According to iNakar et al.l ( 120091 ). in the slow cooling, equation (1T91 leads to Y(j c ,tj) < 
1, then SSC cooling has no effect on the electron distribution. From the analytical solution 
of I Wang et al.l (120101 ). Y{^ GeV , £ 7 ) <Cl can be estimated. With the assumption that the GeV 
luminosity is dominated by synchrotron emission, F synjGeV (13s) ~ 0.83/iJy provide a stronger 
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constraint, -Ek,55 = 0.19A 3 5 S (1 + 2x) 1 / 10 .In such parameter space, we have r 77 <Cl, so the 
absorption does not contribute to the early rising. 

All the above relations show that e ej o, Ek,55 and A 35 , 5 can be expressed with rj 3 and x, 
the physical mechanism for the rising part of the lightcurve can be considered as contribution 
of the structured ejecta. According to observational statistics, the best fitting yields -Ek,55 — 
0.13, A35.5 = 0.02, e efi = 0.3, e B ,o = 5.81 x 10~ 6 , and r] 3 = 0.74. As the result shown in 
Fig. fig:080916C: before the deceleration time, the bulk Lorentz factor of the swept blast 
wave increases with x = 0.5 (which is self-consistent with our assumption for x and thin 
shell case), t cm occurs at ~ 4.5s then electrons are slow cooling, and deceleration begins 
from ~ 6 s, the density jump is expected around 10 5 ' 3 — 10 5 ' 6 s (i.e., 3.25 x 10 18 cm). After 
breaking out to the interstellar medium, the surrounding density is around 1.17 x 10~ 3 cm~ 3 . 



4.2.2. GRB090902B 

Along with the LAT and XRT data mentioned above, our optical photom etry here is 



taken from lPandey et al.l (120101 ). as well as the 8.5GHz data reported by the VLA (Ivan der Horst et al 
20091 ). With the assumption that the late X-ray, optical and radio afterglow and the early 
LAT photons originate from the same region, their spectral indices /3 6 iX = 0.90 ± 0.13 and 
(3ob,o P t = 0.76 ± 0.07 should be consistent with the GeV photons, a b t LAT = 1-4 ± 0.06. The 
same as the above analysis, VQeV > max[i/ c (t/), v m (ti)\ and v c {tn) > ux > v op t > v m (tn), 
respectively, which result in the plausible light curve and spectrum similar to the above case. 
Therefore, roughly speaking, 2.3 < p < 2.6 and 0.17 < e e < 0.43 are the available region. 

Although the peculiar "soft-hard" spectral evolution seems tricky in the prompt emission 
phase, its duration is so short (T — T + 10 s) and it ends before the deceleration time 



( IZhang et al.ll2010l ). which, is understandable in the frame of structured ejecta acceleration 
(we will discuss it in Secj5]), so the external shock model is still applicable. In addition, 
although a 33 GeV photon was detected at 82 s after the trigger, we consider the 11.2 GeV 
photon at ~ 10 s as the highest energy photon of estimation for the lower Lorentz factor 
because the arrival time of the 33GeV photon is far beyond the prompt emission phase and 
may not merely be attributed from the external shock but some other origins. 

The turning point of the LAT light curve occurs at ~7 s. This indicates that a relation- 
ship among E K)55 , r/ 3 , and A 35 5 becomes r] 3 ~ O^i^ggAgg 1 ^! + 2x) 1//4 . Hence, equation 
(ED yields that at t < t dec and f /T 2 = 98 exp(5.94x)(l + x)(l + 2x)- 1 r 3x ~ 2 > 1, so the thin 
shell case is applicable in this fitting. 



Similar to analytical solutions of GRB080916C we mentioned above, characteristic fre- 
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quencies of the photons observed by LAT and later multiband (X-ray, optical, radio) af- 
terglow can be estimated. Additionally, when we come to the radio band, SSA may be 
crucial. 

Here again, we assume that the typical value s(ti) ~ 1/3, while e(tu) ~ 0.0, with model 
parameters (-Ek,55, £b,o, £e,o, and A 35 . 5 ), the characteristic frequencies evolve as 

u c (tj) = 3.29 x 10 9 Hz4^ 5 A 3 - 5 2 5 e^ /2 (l + 2h)^H^\\ + Y( lc ,tj)}- 2 , (20) 

u c {t n ) = 2.71 x 10 9 Hz4^ 5 A 3 - 5 2 5 e- 3 /2 (l + 2x)- 5 / 2 i // 1 / 2 [l + Y{ lc , t n )]- 2 , (21) 

v m {tj) = 6.52 x 10 23 Hzi^e 2 ef Q (l + 2x) 3 / V 8/5 / P 2 , (22) 

y m {tn) = 5.37 x 10 23 Hzi^e 2 e^ 2 (l + 2x) 3 / 2 t / ^ 3/2 / P 2 . (23) 

Setting p = 2.4, e ej o = 0.4 as a trail, the synchrotron emissivity is given by 

i^n,Gev(10s) = 14.58/i JyE^e^ (1 + 2x) 3 / 10 [l + Y( 7c , lOs)]" 1 . (24) 

F syntX (1.17 x 10 5 s) = 6.72 x 10 5 fiJjE 1 ^ A 35 . 5 e^ 20 '(1 + 2x) 31 / 20 ~ 0.18/iJy (25) 

From equation ([25]), we have e B ,o = 1.87 x lO- 8 ^ 1 ^" 2 ^ 17 ^ + 2x)~ 9 / 5 , so t cm < 10s in the 
case that A 35 , 5 > 10~ 2 . Therefore, i/ m (10s) < z/ c (10s) and F syni G e v(10s) < 1.59/iJy, which 
yields, 

1.54E K ,55^35 3 f (1 + 2x) 3 / 25 < 1 + F( 7c , 10s) < 55E k ,5 5 aS(1 + 2x)- 9 / 50 . (26) 

Meanwhile, for the radio band at tin s, because of the light curve index ~ 1/3, SSA can be 
neglected, we have 

z/ m (5.23 x 10 5 s) = 9.14 x 10 10 A^Jj| /17 (l + 2x) 10/17 > 8.5GHz, (27) 
^n,8.5GHz(5.23 x 10 5 s) = 1.89 x 10 2 /xJyA 3 / 5 5 (l + 2x)~ 3 / 5 ~ 56^/iJy. (28) 

Therefore, it should be 0.05 < A 35 . 5 (l + 2x)^ < 0.24, and if A 35 . 5 (l + 2x)~ 1 t I < 46.7 
(e.g.A 35 . 5 (l + 2x)~ x < 0.24 at 10s or 0.05 < A 35 . 5 (l + 2x)^ ~ 0.05 at 10 3 s), 

= 1-8 x 10^ K i 5 A 2 / 5 5 (l + 2x)-/ 5 t 3 / 5 [l + F( 7c ,t / )] < 0.10A 3 / 5 5 (l + 2x)- 3 /V /5 < 1. 

(29) 

Under such a parameter space, equation (|2"9"|) leads to Y c tl < 1, SSC makes a little 
effect because of the Klein-Nishina suppression, so the estimation that Y{ r )GeV)^i) is 
self-consistent. With the assumption that the GeV luminosity is dominated by synchrotron 
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emission, F S j / „ i Gev(10s) ~ 1.59/iJy provide a stronger constraint, -Ek,55 = 0.65A 3 5 5 (1 + 
2x)~ 2 / 17 and 

Using the above relations, we can express e e ,o, Ek,55 and A 35 . 5 by 773 and x, and fit the 
data according to observational statistic. The best fitting requires -Ek,55 = 0.44, A35.5 = 0.06, 
e 6i o = 0.4, eBfl = 4.56 x 10~ 7 , and r] 3 = 0.60. Fig. fig:parameter shows the circumburst 
density, radius and characteristic frequency, and Compton Y Ge y and Y lc factor evolution 
under such a parameter space. Fig. fig:090902B shows the spectrum and light curve fitting. 
Before the deceleration time (7s), the bulk Lorentz factor of the blast wave increases with a 
slope of x = 0.8 (which is self-consistent with x and thin shell case). The shock-accelerated 
electrons are cooling fast early before ~ 0.72 s, the density jump is estimated around 10 5 ' 1 — 
10 6 s (i.e., at radius of 4.48 x 10 18 cm). After breaking out to the interstellar medium, the 
surrounding density is as low as 1.91 x 10~ 3 cm -3 . These parameters are consistent with 
the radio light curve except for the first data and is reliable when the jet break takes place 
nearly at the same time. We find that the jet angle is 0.04 rad and collimation-corrected 
energy is Ej et ~ 1.51 x 10 51 ergs. 



4.2.3. GRB 090926 A 



V band and R band data of this burst are taken from ISwenson et al.1 (120101 ) and 



Rau et al.l ( 120101 ). respectively. Our analysis procedure is the same as mentioned above. 
Due to the hypothesis of the same source, the temporal index of the LAT photons oiob^LAT = 
2.05 ± 0.14, should be consistent with the spectral indices (3 b )X = 1.12 ± 0.13 and (5 b, op t = 
1.03±0.05. This leads to a fairly broad plausible range for 3 > p > 2 and e e ~ 10 _1 (e.g. under 
the condition of VQeV > rnax[z/ m (t/), u c (tj)], v c (tn) > v^y requires for 3 > p > 2.96 and 
0.62 > e > 0.37, while 2.5 > p > 2 and 0.70 > e > 0.45 as long as v c {t n ) < v x )- However, 
since the light curve after the density jump is a little bit steeper than the early part, e (the 
latter the smaller) should not deviate too much from the typical value (i.e., 1/3, otherwise, 
the decrease of loss efficiency will cause the decaying slope much shallower), whereas the 
decay of LAT is steep, which requires a large p (i.e., 3.0). Besides, in the afterglow phase, 
we can reasonably believe that the first flare at 70 — 95 ks comes from an interaction of the 
blast wave and density jump. The second "flare" may come from late energy injection or 
another density jump, and in this paper, we smooth it as the simple power law. 

Because the deceleration time is estimated as ~ 16 s, 773 can be written by i] 3 ~ 

0. 27 J E^ 5 A 3 " 5 1 / 4 (l+2x) 1 / 4 . Hence, att < t dec , f/T 2 = 512 exp(8.32x)(l+x)(l+2x)^ 1 r M > 

1, so the reverse shock is negligible in this fitting. With the typical value e(tj) ~ 1/3, while 



- 19 - 



£{tn) ~ 0.0, The characteristic frequencies evolve as 

i/ e (f,) = 3.09 x 10 9 Hz4^ 5 A 3 - 5 2 5 e B 3 /2 (l + 2h)-^H i 2 / 5 [1 + Y( lc ,tj)]- 2 , (30) 

u c (tjj) = 2.34 x l0 9 RzE^ 2 55 A- 5 2 5 e~f(l + 2x)- 5 / 2 t 7 // 2 [l + Y( lc ,t n )}- 2 , (31) 

u m (t z ) = 7 A3 x 10 23 Hz^e 2 e^ /2 (l + 2x) 3 / V 8/5 / P 2 , (32) 

v m (t u ) = 5.63 x 10 23 HzEjg 5< £ e^(l + 2x) 3 / 2 t // ~ 3/2 / P 2 . (33) 

Setting p = 3.0, e ei o = 0.4 as a trail, the synchrotron emissivity is calculated by 

iV,Gev(36s) = 28.7 /U JyE^ 4 55 e^(l + 2x) 3 / 4 [l + Y( 7c , 36s)]- 1 , (34) 

F syn , x (5.3 x 10 4 s) = 1.42 x 10 6 /iJyEK,55A 3 5.5eB,o(l + 2x) 2 ~ 0.33/xJy. (35) 

From equation (|55]h we have e B ,o = 2.32 x lO -7 ^ 1 ^^^ + 2x)^ 2 , so t cm = 4.31/£ K ,55 < 
36s in the case that E K>55 > 0.12. Therefore, z/ m (36s) < z/ c (36s) and -F S j /n ,Gev(36s) < 0.37/zJy, 
which yields, 

l.71E Kj55 A;l f 5 \l + 2X) 1 / 4 < 1 + F( 7c , 36s) < 8.38E K , 55 . (36) 

Hence,as long as A^ 2 5 (l + 2x)- 1 / 2 t / 3 / 5 < 65.8 (e.g., A 35 . 5 (l + 2x)-* < 58.7 at 36s or 
A 35 . 5 (l + 2x)- x < 1 at 10 3 s), 

= 1-8 xlO- 3 ^^ < 0.015A 3 / 2 5 (l + 2x)- 1 /V/ 5 < 1. 

lc{tl) 

(37) 

In the slow cooling case, equation f|3T|) leads to Y(^ c , ij) < 1, SSC makes a little effect 
due to the Klein-Nishina suppression, because the estimation of Y{pf GeVl tj) <Cl is self- 
consistent. 

Meanwhile, during 1 < tj < 10 2 ' 5 s, in the case of fast cooling (£55 < 0.12) that 
z/ m (36s) > z/ c (36s) and F syn:GeY (3Qs) < 0.37/iJy, A 35 . 5 (l + 2x) x /s < .026, 

— O QS \y m3/i-i/2 



lm{tl) 



2.38 x lO 3 ^ 2 (1 + 2x) 1/2 tr 8/5 > 1, (38) 



According to iNakar et al.l (120091 ). this is within the strong KN regime (case II). Analytical 
discussion with the above confines yield to the possibilities of case lib and lie, both leads to 
Y (jGeViti) «d. Hence, the contribution from SSC could be neglected, too. In both fast and 
slow cooling, i r S yn i Gev(36s) ~ 0.37/uJy provide a stronger constraint, E K 55 = 0.70A35 4 5 (1 + 
2X)" 1 / 4 and r 77 < 1. 
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Using the above relations, we can express e e> o, -Ex,55 and A 35 , 5 by 773 and x, and fit 
the data according to observational statistic. When E K 55 = 0.43, A35.5 = 0.03,e ej0 = 0.35, 
e B,o = 8.18 x 10 -6 , and 773 = 0.56, fitting seems plausible (Fig. fig:090926A). Before the 
deceleration time, the Lorentz factor of the blast wave increases with a slope of x = 0.3, 
tcm ~ 8.8 s and from ~ 16s significant deceleration appears. The blast wave enters the ISM 
around 10 4 ' 8 — 10 4,9 s (i.e., at radius of 4.41 x 10 18 cm). The medium density surrounding 
the bubble is 1.3 x 10~ 3 cm~ 3 . 



4.3. Analysis of Results 



Tab. [2] shows all the best fit parameters for the above three GRBs. Common results can 
be found: 

(1) In all the three burst fittings, the thin shell case is applicable, the very beginning 
of GeV flux densities are forward-shock dominated. The rising slope comes from the con- 
tribution of the power-law increasing bulk Lorentz factor of structured ejecta. These are 
self-consistent with the thin shell case, indicating that the central engine is very likely to be 
the core collapse of a massive star (we will discuss it in Secj5]). Our model implies that the 
central object would keep on accreting matter even after the LAT trigger. 



all 



(2) The remaining isotropic kinetic energy of the fireball after the prompt emission are 
- 10 54 ergs, indicating ~ 30% — 90% of the total energy has been emitted in the prompt 
emission. This is because t cm occurs before the deceleration time for all the three GRBs 
and the following dynamics could be considered as being quasi-adiabatic. As a matter of 
fact, e e in our sample fittings is ^ 2/ 3, which meets t he requirement of radiation efficiency 
related to the analytic procedure from lWu et al.l (120051 ). Practically, slow cooling of electrons 
starts before the deceleration time, and a jet break appears much later, so Wu's spherical-like 
solution is tenable in this piece. A low value of (~ 10~ 7 — 10 -5 ) indicates a weak magnetic 
field, especially for GRB 80916 C and GRB 090902B. These GRBs have been studied by 
Kumar fc Barniol Duran I (120091 ). whose inferred magnetic field (along with several other 
LAT-detected events) is consistent with shock compression of a modest circumstellar field 
(B £3 30/uG). That is to say, no dynamo process is necessary to generate the magnetic 
field needed for the observed synchrotron afterglow emission. In our best fittings, we find 
B > 380/iG for GRB080916C and B > 106/xG for GRB09 0902B ev en wh en T + 10 7 s, being 
broadly consistent with the results of lPiran &: Nakarl (120101 ) and lLi I (120101 ). Meanwhile, when 
the blast wave enters the const ant- density circumburst medium, the preshock magnetic field 
can be calculated with our fitting parameters as B = (27rm p e,Bn) 1//2 ~ 8.0/iG for 080916C 
and 2.6/iG for 090902B, that suggests no magnetic field amplification. 
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(3) Although the ratio of e e and e# is large, synchrotron self-Compton is not significantly 
contributed to the observed flux density in the GeV band because of the Klein-Nishina 
suppression effect. In our sample fittings, only synchrotron radiation from the external 
forward shock is enough to contribute to the high energy emission and the synchrotron self- 
Compton component can be neglected because Y{^fGev) "C 1 in this piece. In particular, for 
GRB 090902B, from our best-fitting parameters, the maximum electron Lorentz factor at 
82s indicates the highest photon energy can be as high as ~ 37 GeV, so it is plausible that 
33 GeV photon may result from synchrotron radiation. 

(4) The wind parameter for thre e longest bursts is A35.5 ~ 10~ 2 — 10 _1 , which is reliably 
low. According to bai fc wJ fl2003h . the progenitor star of GRB 080916C, GRB 090902B, 
and GRB 090926A could have been in a cloud. Due to a high pressure of the cloud, a 
speedy wind will be slowed down by a pair of shocks (viz, a reverse shock that propagates 
into the wind gas and a forward shock that propagates into the cloud). This interaction 
produces a stellar wind bubble and forms a density-jump at a radius of ~ 10 18 cm. The 
homogeneous density outside this jump is estimated to be in the range of n ~ 10~ 3 cm~ 3 , 
which is much less than the typical density of an interstellar medium, but could not be 
excluded from the possibil ity of existence in areas of active star formation as the interiors of 
a prexisting superbubble (IPanaitescu fc Kumarll200ll ; IScalo &: Wheeled l200ll ). If confirmed 
by the superbubble observation, this can provide a piece of evidence for a connection between 
long-duration GRBs and broad-lined SNe Ib/c. and would be natural to expect that massive 
progen itors of GRBs explode where we obser ve star formation. Meanwhile, in other fitting 
works (jCenko et all 120101 ; iLiu fc Wand |2010| ) , they have also confronted with the problem 
of low density, and suggest the lower metallic ity progenitors with minimal pre-explosion 
mass loss or selection effects (jCenko et al.ll2010f). Hence, alt hough the LAT emission at least 
during T 90 , from the spectral perspective (jZhang et al.ll2010l ). is likely to connect to the GBM 
emission of internal origin, the fitting facts tentatively suggest a possibility that the LAT 
observation records an internal-external shock transition and at least in the LAT decaying 
phase the external forward shock emission is dominated. 



5. Conclusions and discussions 

It is widely believed that long duration gamma-ray bursts (GRBs), like hydrogen- 
deficient Type Ib/c supernovae (SNe Ib/c), result from the core collapse of a massive star. 
The main characteristic of setting GRBs apart from other SNe is that a substantial frac- 
tion of the explosive energy is coupled to relativistic ejecta. A compact central engine is 
responsible for accelerating and collimating a jet-like outflow and driving a SN explosion 
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( IWoosley fc Blooml 120061 ; iGehrels et al.l 120091 ; ISoderberg et al.l |2010[ ) . The precise nature of 



the central engine which powers GRB-SNe, however, remains an open question. 

In this paper, we have undertaken extensive broadband (>100 MeV, radio, optical, 
and X-ray) observations of three long-duration GRBs (GRB 080916C, GRB 090902B, GRB 
090926A) detected by the LAT instrument on the Fermi satellite. The bulk Lorentz factors 
imply that the fireball is ultra-relativistic, which indicates that hyper-energetic bursts carry 
as high as E iso ~ 10 54 ergs in the blast wave. The temporal indices a and spectral indices 
(3 in Tab. [TJ which are consistent with the observed simple power law slopes, show that 
our detailed density transition consideration can reconcile the plateau and shallower decay 
in the later afterglow. Tab. |2] displays the estimated physical parameter space of the best 
fit, which is all reasonably certificated by the theory mentioned in SecJU Specifically, the 
wind parameter A 35 5 ~ 10~ 2 — lO^ 1 s hed us light on the central engine of stellar collapsars 



and their associatio n with SNe Ib/c (IMacFadyen fc Woosleylll999l ; IWoosley fc Hegerl 12006 



Modiaz et al.ll2008T ). 



Furthermore, the tricky yet popular early rise of the LAT (>100MeV) light curve 
may indicate whether high energy photons share the same source as the low energy af- 
terglows. Aside from two main kinds of explanation we mentioned above, leptonic or 
baryonic models, an alterna tive possibility comes from the high-latitude prompt emission 
( Kumar &: Panaitescu 200ol ) and thus the flux density evolves as F v {t) oc t _2_/3 . However, 



according to P> b,LAT in Tab. [TJ, the calculated temporal index is too steep to explain the early 
high energy emission decay. Or else if the GeV component results from the same conven- 
tional unstructured ejecta as the late X-ray and optical afterglows, before the deceleration 
time (t < tdec, that is, we only consider the thin shell case for long bursts), in the early 
wind bubble n oc R~ 2 , the bulk Lorentz factor is nearly constant, and thus the electronic 
cooling and minimum frequencies evolve as v c oc t and z/ m ; n oc t~ l , the peak flux density 
F max oc t°. Even when the SSC makes a significant effect, the IC electronic cooling and 
minimum frequencies and IC peak flux density evolve as fic,c oc t 3 and fic,m oc t^ 1 and 
F/c,max oc t~ l . Since max(z//c,c, ^ic,m) > ^Gev > max(i/ c , z/ m ) around the peak of LAT light 
curve, F syn ^ GeV oc £( 2 ~p)/ 2 (for 3 > p > 2, this slope is not enough), as well as F ICy GeV oc i~ 2 / 3 
or t~(p +1 )/ 3 for slow cooling and Fjc,GeV oc t~ 2 or t~ 5//2 for fast cooling. No matter whether 
take the KN cutoff into account, this might be inconsistent with the observed rise of the 
LAT light curve. 

Our explanation for the initial rising part is simply based on the assumption that T oc 
t H . Therefore, v c oc t* +x , v m oc and F syn , GeV oc e -^it {1+ ^ 2 ~ p)/2+2x when v GeV > 

max(i/ m , u c ). Co-contribution of the optical depth and structured ejecta bring about the 
rising. This assumption can be understandable within the framework of of the core collapse 
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of a massive star (i.e., SNe Ib/c). If the central engine ejects matter continuously, a clean 
tunnel would be left after the early ejecta sweeps up the stellar envelope, and the succeeding 
ejecta, whose energy is attained via continuous accretion even after the LAT trigger and 
whose Lorentz factor is larger, catches up with the early ejecta and injects more energy into 
the blast wave, giving rise to a rise of the LAT light curve. On the other hand, the thermal 
component can also be explained under such assumption that the previous jet lost kinetic 
energy due to impedance from sho ck frontier, then turned it into thermal radiation. This 
hypothesis have been provided by (IMacFadyen et al.l 120011 ) that, in the case of red-giant or 
blue super giants , the powerful jet lose its energy input at its base in the course of overtaking 
the sub-relativistic weak supernova shock before reaching the surface, and eject little highly 
relativistic matter; While in a helium star, the jet breaking out with continuously accretion 
power receiving at t he base will lead its motion become highly relativistic . La t er, so me other 
possibilities, such as IZhang et al.l (120031 ) 's baryon pollution scenario and iFanl ()2009|)'s initial 
envelope choked model, also result in the earlier power-law rising distribution of bulk Lorentz 
factor. 

In our frame, sub-MeV photons provide the background annihilation with the later 
GeV ones. On one hand, photon with the largest energy f r om th e external shock fasten the 
constraint of lower Lorentz factor of the ejecta [in lGranot I (120101 ). the lower limit of the bulk 
Lorentz factor should be > 1000]; on the other hand, even the ejecta is slower than the lower 
limit, it is possible that the earlier rising is paretically contributed from opacity effect. 

About the observed time delay of the pulses and variabilities in both LAT and GBM 
bands, some other modulation effects, except for optical depth can account for the concur- 
rent variabilities in both bands, such as inhomogeneous surroundings, energy in jection, etc.. 



Thes e factors indicate that the variability connection between LAT and GBM (IZhang et al. 



20101 ) cannot give a sufficient clue to the internal/external origin of the GeV photons. 



In particular, as for GRB090902B, one of the tricky properties is its soft-hard-soft quasi- 
thermal spectrum evolution, since the soft-hard evolution appears within T90, it can be better 
understood by our model that the initial outflow was dissipated as thermal component by the 
denser shock front/photosphere, when it breaks out of the envelope, the consequent highly 
relativistic ejecta in a clean tunnel can be dissipated strong enough to produce energetic 
non-thermal emission, even the photosphere can be outshined to some extent, therefore, the 
spectrum evolves from soft to hard. 



This work is supported by the National Natural Science Foundation of China (grants 
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Table 1: According to F u oc t~ a v~", cx^at, &x,o,i> a x,o,2 are the estimated temporal indices of the LAT band photons, 
X-ray ± optical photons in the wind bubble, and X-ray + optical photons in the homogeneous medium from our best-fit 
parameters. a b t LAT, &ob,x, &ob,o are the observed temporal indices of GeV photons, simple power-law X-ray photons and 
optical photons. Plat and /3 X>0 are the estimated spectral indices of the LAT band photons and X-ray + optical photons 
in either medium from our best-fit parameters. (3 a b,LAT, (3 b,x, and (3 b >0 are the observed spectral indices of G eV photons, 
simple power- law X-ray phot ons and optical photons. High energy data are taken from lZhang et al.l ( 120101 ) . X-ray data 
are taken from Swif t Group (jEvans et al.l 120071. 120091) . In addition , the data labeled with su perscript a are taken from 
Greiner et al.l (120091 ). the data with b are from lPandey et al.l (120101 ). the data with c are fr o mlSwenson et al.l (120101 ). — * 



represents no data available. For 090926A, the data labeled with * are from lCenko et al.l (120101 ). the temporal slope of 
the X-ray and optical band in iRau et al.l (120101 ) is: a = 1.6 ± 0.2 before the first flare, a = 1.75 ± 0.04 after the second 
flare, and a = 1.63 ± 0.01 between, which is consistent with our estimated value. 
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Table 2: The left part of the table are best-fit parameters for GRB080916C, GRB090902B, and GRB090926A, the right 
part are prompt emission energy fraction, radius when a density jump occurs and circumburst density after the jump, 
respectively. 
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Fig. 1.— Best fit for GRB 080916C: p = 2.32, x = 0.5, Egfc = 0.13, r] 3 = 0.74, A 35 , 5 = 0.02, 
e ej o = 0.3, €b,o = 5.81 x 10~ 6 , e = 2.1, tfe rea fc = 10 5 ' 3 s, t en d = 10 5 ' 6 s, no jet break appears 
during the modeling period 
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Fig. 2.— Best fit for GRB 090902B: p = 2.4, x = 0.8, E%° 55 = 0.44, rj 3 = 0.60, A 35 . 5 = 0.06, 
e e , = 0.4, e Bi0 = 4.56 x lO" 7 , e = 2.0, t break = 10 51 s, t end = 10 6 s, t jet = 10 6 s 
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Fig. 3.— Best fit for GRB 090926A: p = 3.0, x = 0.3, Egfe = 0.43, rj 3 = 0.56, A 35 . 5 = 0.03, 
e 6i0 = 0.35, €b,o = 8.18 x 10~ 6 , e = 2.5, t break = 10 4 8 s, t en d = 10 4 9 s, no jet break confirmed 
until 21 days 
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Fig. 4. — Physical parameter evolution with the best fit parameter for GRB090902B, similar 
to the other two bursts. Upper-left panel: a. circumburst density with radius, upper-right 
panel: b. blast wave bulk Lorentz factor with time; lower-left panel: c. Compton Y factor 
(<Cl due to KN effects) with time; lower-right panel: d. characteristic frequency evolution. 



